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Abstract 

We study integrable realizations of conformal twisted boundary conditions for si{2) 
unitary minimal models on a torus. These conformal field theories are realized as the 
continuum scaling limit of critical G = A, D, E lattice models with positive spectral 
parameter u > and Coxeter number g. Integrable seams are constructed by fusing 
blocks of elementary local face weights. The usual A-type fusions are labelled by the 
Kac labels (r, s) and are associated with the Verlinde fusion algebra. We introduce 
a new type of fusion in the two braid limits u — > ±ioo associated with the graph 
fusion algebra, and labelled by nodes a,b £ G respectively. When combined with 
automorphisms, they lead to general integrable seams labelled by x = (r, a, b, k) £ 
(v4 9 _2, H, H, Z2) where H is the graph G itself for Type I theories and its parent 
for Type II theories. Identifying our construction labels with the conformal labels of 
Petkova and Zuber, we find that the integrable seams are in one-to-one correspondence 
with the conformal seams. The distinct seams are thus associated with the nodes of the 
Ocneanu quantum graph. The quantum symmetries and twisted partition functions 
are checked numerically for \G\ < 6. We also show, in the case of D21, that the non- 
commutativity of the Ocneanu algebra of seams arises because the automorphisms do 
not commute with the fusions. 



1 Introduction 

There has been much recent progress [l|, IE IE 0, IE IE S El; on understanding integrable 
boundaries in statistical mechanics, conformal boundary conditions in rational conformal 
field theories and the intimate relations between them on both the cylinder and the torus. 
Indeed it appears that, for certain classes of theories, all of the conformal boundary condi- 
tions on a cylinder can be realized as the continuum scaling limit of integrable boundary 
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conditions for the associated integrable lattice models. For s£(2) minimal theories, a com- 
plete classification has been given of the conformal boundary conditions on a cylinder. 
These are labelled by nodes (r, a) of a tensor product graph A £g> G where the pair of graphs 
(A, G), with G of A-D-E type, coincide precisely with the pairs in the A-D-E classification 
of Cappelli, Itzykson and Zuber 0| . Moreover, the physical content of the boundary condi- 
tions on the cylinder has been ascertained |4|, |lCj by studying the related integrable boundary 
conditions of the associated A-D-E lattice models llj for both positive and negative spectral 
parameters, corresponding to unitary minimal theories and parafermionic theories respec- 
tively. Recently, the lattice realization of integrable and conformal boundary conditions for 
N — 1 superconformal theories, which correspond_to the fused A lattice models with positive 
spectral parameter, has also been understood 
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In this article, we use fusion to complete the task of constructing integrable real- 
izations of conformal twisted boundary conditions on the torus. Although the methods are 
very general we consider s£(2) unitary minimal models for concreteness. The key idea is that 
fused blocks of elementary face weights on the lattice play the role of the local operators 
in the theory. The integrable and conformal boundary conditions on the cylinder are con- 
structed [4] by acting with these fused blocks on the simple integrable boundary condition 
representing the vacuum. For the usual A-type fusion, associated with the Verlinde fusion 
algebra, this leads to integrable seams labelled by the Kac labels (r,s) G (A g ^ 2 , A g ^\). In 
this paper we introduce a new type of fusion of G-type related to the graph fusion algebra. 
Integrable seams of this type are labelled by (r, a) G (A g _ 2 , H) where H is the graph G itself 
for Type I theories and its parent for Type II theories. By the generalized Yang-Baxter 
equations, these fused blocks or seams can be propagated into the bulk without changing 
the spectrum of the theory. The seams so constructed provide integrable and conformal 
boundary conditions on the torus. Fixed boundary conditions a G G on the edge of the 
cylinder are propagated into the bulk by the action of the seam (l,a) on the distinguished 
(vacuum) node 1 G G. Lastly, automorphism seams, which play no role on the cylinder, play 
a crucial role on the torus by providing the extra label giving rise to the complement of the 
left and right chiral subalgebras in the Ocneanu graph. 

In general, for rational conformal field theories on the torus, we expect the two types of 
fusions supplemented by the automorphisms to generate all of the integrable and conformal 
seams. In this paper we discuss this assertion in the context of the A-D-E unitary minimal 
models. 

The paper is organized as follows. In Section |2] we define the A-D-E series, giving their 
graphs (Sec. l2.ljh their (proper or improper) graph fusion algebras (Sec. 12.2)1 . their Ocneanu 
graphs (Sec. 12. 3|) and their associated twisted partition functions (Sec. 12.31 and l2.4|) . The 
presentation is self contained. In Section|3]we describe the lattice realization of these twisted 
boundary conditions. In particular, we define the A-D-E lattice models (Sec. l3.1Jl . the fusion 
projectors (Sec. I3.2p . the associated fused faces (Sec. I3.3p and the integrable seams (Sec. 3.4 
to 3.6). We construct the transfer matrices in Section 0] composed of regular faces and seams. 
This is described for the single-row transfer matrix on the torus (Sec. 4.1) and the double-row 
transfer matrix on the cylinder (Sec. 14.3(1 . The spectra of the transfer matrices and finite-size 
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corrections are described in Section |3J The free energies are computed (Sec. I5.2J) and the 
numerical conformal parts are identified with the twisted partition functions (Sec. 15.31 to 
5.6). 

2 A-D-E Fusion Graphs and Partition Functions 

A-D-E classifications appear in a variety of contexts, namely, graphs, solvable lattice models, 
s£(2)k (Wess-Zumino-Witten) models at level k, and s£{2) minimal models. 

2.1 A-D-E Graphs 

The basic A-D-E objects are graphs. A simple graph G is given by its vertices (or nodes) 
a G Go and edges (a, b) G G\ C Go x Go- We are concerned with unoriented graphs, 
(a, b) G G\ => (b,a) G G\. The A-D-E graphs, which are the Dynkin diagrams of simply 
laced Lie algebras, are presented in Table H The number g is the Coxeter number of the 
graph G and the exponents Exp(G) are a subset (with multiplicities) of the nodes of the A L 
graph sharing the same Coxeter number as G. 

Graph G g Exp(G) Type/if T 

A L I I \ - L . L+l 1,2,-.. I Z 2 

<i+i 
2£+2 1,3,-.. ,2£+lJ + l I Z 2 

e+2 

<e+i 
21 + 2 1,3,... ,2£+l,£ + l ll/A 2e+1 Z 2 

1+2 

6 

E 6 llli 5 . 12 1,4,5,7,8,11 I Z 2 

3 

7 

E? i — S— 2 — I — S — S 18 1,5,7,9,11,13,17 \\/D m 1 

8 

£ s I I I j I g I 30 1,7,11,13,17,19,23,39 I 1 



Figure 1: A-D-E graphs corresponding to the Dynkin diagrams of the classical 
A-D-E simply laced Lie algebras. The nodes associated with the identity and the fun- 
damental are shown by *, □ respectively. Also shown are the Coxeter numbers g, exponents 
Exp(G), the Type I or II, the parent graphs H ^ G and the diagram automorphism group 
T. The D 4 graph is exceptional having the automorphism group S3. 

A graph G is completely encoded by its adjacency matrix which we denote by the same 
letter G. It is a symmetric non- negative integer square matrix whose rows and columns 
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are labelled by the vertices of G, defined by G a b = 1 if a and b are adjacent and G a b = 
otherwise. What is so special about the A-D-E graphs is that (along with the tadpole 4 series) 
they are the only ones whose spectra lies in the open interval (—2, 2). The Perron- Frobenius 
theorem implies that the largest eigenvalue of these adjacency matrices is non-degenerate, 
real and positive and its eigenvector can be chosen to have non-negative entries. They are 
given explicitly in terms of g-deformed integers S n = [n] q = ^qpi with q = exp(ni/g). The 
largest eigenvalue is S 2 = [2] g and the eigenvector t/j is 



Kk<L 



*P E6 = ([l] g , [2] q , [3} q , [2} q , [l] q , (2.1) 

^ = ([l] g , [2] g , [3] g , [4]„ |k, |k, 

^ = ([1] ? J2] 9 J3] ? J4] 9 J5] ? ,||,||,||). 

2.2 Graph fusion algebras 

The integer linear span of the nodes of the graph can be given a structure of a commutative 
graph fusion algebra. We first specify two vertices, the identity vertex * and the fundamental 
vertex □. They are indicated in Table Q and are respectively the vertices labelled 1 and 2 
in the Al, D 2 £+2, Eq and E 8 cases, known as Type I theories, the vertices 21 + 1 and 21 — 1 
for the D21+1 case, and the vertices 1, 2 for the E7 case, known as Type II theories. The 
former models give rise to proper graph fusion algebras with non negative integer structure 
constants or Non-negative Integer Matrix Irreducible Representations (nimreps) which are 
understood as adjacency matrices. The latter have some negative structure constants and 
do not form proper graph fusion algebras. 

The algebra is defined by stating that the edges of the graph G encode the action of the 
fundamental element □: 

a □ =J2 b - ( 2 - 2 ) 

The identity gives one row of the algebra table, the previous formula gives another, commu- 
tativity and associativity determine the rest. On the -D4 example, □ = 2 and associativity 
gives 

(4-3)(D □) = (4 D-3 □)□ = () 

= (4-3)(l + 3 + 4) (2.3) 

r 3 3 = 4 

so that 4 — 3 = 33 — 44 and □ 3 3 = □ shows that < and the expansion of 3 3 3 



implies 3 4 = 1. 



44 = 3 



4 The tadpole graph Tl is obtained from the graph Al by adding a loop at the final vertex; it is not a 
simple graph. 
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The structure constants of this algebra are denoted N 

ab = Y,N ab c c. (2.4) 

ceG 

The definition of the algebra implies iV* = I and Nu = G and these matrices themselves 
form the regular representation of the algebra with the usual matrix product 

N a N b = Y,N ab c N c . (2.5) 

ceG 

As it is a commutative algebra containing the adjacency matrix, its common set of eigen- 
vectors is given by an orthogonal basis of eigenvectors of G. They are labelled by Coxeter 
exponents (we have only given the Perron- Frobenius eigenvector) and the spectral decompo- 
sition of each matrix onto its eigenvectors give these integers through a Verlinde-like formula 

jeExp(G) V* 

Some algebra tables are given in Tables HHS1 In the case of the graph Al this reduces to 
the usual Verlinde formula, the structure constants are denoted by Nij k and the matrix of 
eigenvectors by S. For an A-D-E graph G with Coxeter number g — L+ 1, another algebra 
of non-negative integer matrices with the Al structure constants Nij k is given by the fused 
adjacency matrices ni defined by the s£{2) recurrence relation 

nx — I, n 2 = G, n i+1 = n 2 ni - Hi_i for 2 < i < g — 1, (2.7) 

which closes with n, = for i > g — 2 and 

_ f I, for D 2i , E 7 , E 8 
% " 2 " 1 a, for A L , D 2 ^ E 6 ^ 

where a is the Z 2 graph automorphism. Clearly, E-j and E$ do not admit a Z 2 automorphism. 
So the fusions contain the Z 2 graph automorphism in all cases where it exists except for D 2 £. 
The matrices rii also satisfy a Verlinde like property 

jSExp(G) b * 

with the algebra structure 

mnj = N H kn k, riiN a = ^n ia b N b . (2.10) 

keA L beG 

The matrices are in fact linear combinations of the graph algebra matrices -/V a 

ni = J2 n n a N a . (2.11) 

aeG 
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Zi = 


XiXl 


+X2X*2 


+X3X3 


+X4XI, 


2 


2 


1 + 3 


2 + 4 


3 


Z 2 = 


X2X*1 


+ (Xl+X3)X*2 


+ (X2 + Xl)X*3 


+X3Xt 


3 


3 


2 + 4 


1 + 3 


2 




X3X*1 


+ (X2 + Xi)X*2 


+ (Xl + Xs)Xs 


+X2Xi, 


4 


4 


3 


2 


1 


Z 4 = 


XaX\ 


+X3X*2 


+X2X*3 


+XlX*4- 



Table 1: Fusion matrices, graph fusion algebra and twisted partition functions of A4 in terms 
of afflne s£(2) characters Xs- 
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0100 
0010 
0001 



No = 



0100 
1011 
0100 
0100 



N» = 



0010 
0100 
0001 
1000 



n = 



Zy — Z 3 i — Z41 — X2X2 

Xl = Xl + X5, X2 = X2 + X4, X3 = XA = X3 
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4' 
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2' 
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1' 
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1 
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10 
10 
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Z l = 


X1X1 


+X3X3 


+X4X4, 


2 


2 


1+3+4 


2 


2 






+X2X3 


+X2X4 — z\ m -, 


3 


3 


2 


4 


1 


Zz = 


XsXt 


+X4X3 


+X1X4 = ^4, 


4 


4 


2 


1 


3 


Z 4 = 


X4X*1 


+X1X3 


+ XsXl 




£>4 

Table 2: Fusion matrices, graph fusion algebra, twisted partition functions, Ocneanu algebra 
and Ocneanu graph of D 4 . The extended chiral and ambichiral subalgebras are bold. 



N 1 = 




N 2 



-10 10 0' 
1 1 
10 10 
10 
1 0, 



N.= 



10 
10 10 
10 11 
10 1 
110 



N 





Table 3: The D 5 graph fusion algebra is not a proper graph algebra. 
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iVi 



1 0- 
10 
10 
10 
10 
■0 1 



No. 



1 0- 
10 10 
10 10 
10 11 
10 
1 0- 



iv 3 



■0 1 0- 
10 10 
10 10 11 
10 2 
10 1 

•001010- 



Na 



■0 1 0- 
10 11 

10 2 

1 2 1 1 / J 
10 10 

•010100- 



Ne, 



1 0- 
10 
10 1 
10 10 I ' 
10 10 
1 0- 



N R 



1 
10 
10 10 
10 10 
10 
10 1 





1 


2 


3 


4 


5 


6 


1 


1 


2 


3 


4 


5 


6 


2 


2 


1+3 


2+4 


3 + 5 + 6 


4 


4 


3 


3 


2+4 


1+3+5 + 6 


2+4+4 


3 + 6 


3 + 5 


4 


4 


3 + 5 + 6 


1+3+3+5+6 


2+4 


2+4 




5 


5 


4 


3 + 6 


2+4 


1 + 5 


3 


6 


6 


4 


3 + 5 


2+4 


3 


1 + 6 



Z\ = XiX*i + X3X3 + X5X5 + XeXt, 

Z2®\ = X2XI + (X 2 + X 4 )X3 + X 4 (Xb + Xe)* = Z{ m , 

Zz = XsXl + (Xi + X3 + X 5 + X 6 ) Xl + (X 3 + Xs) Xl + (X 3 + X 5 ) Xl, 

Za®\ = X4X1 + (x 2 + 2 X 4 ) Xl + (x 2 + X 4 ) (X5 + Xe)* = ^i® 4 > 

^5 = X 5 xt + (X 3 + Xe) Xl + (Xi + Xs) X5 + X 3 Xe\ 

^6 = X6Xl + (X3 + X 5 )X3 + X3(Xl+X 6 )X5 + (Xl+X6)X6> 

^1' = X2X2 + X 4 X 4 , 

Zv = lx 2 + x 4 | 2 + IxJ 2 , 

Zy = Z & = x 2 xl + X4X2 + lx 4 | 2 - 

Xl = Xl + X9, X2=X2 + X8, X3 = X3 + X7, 
Xi = Xi + X6, X5 = X6 = X5- 

Table 4: The graph fusion algebra of D$ and its twisted partition functions in terms of 
extended characters. The extended chiral subalgebra T is shown bold (see Fig. 0J) and the 
extended characters are given in terms of the A$ characters. 
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1 0- 
10 
10 
10 
10 
•0 1 



No 



1 0- 
10 10 
10 10 1 
10 10 
10 
1 0- 



1 0' 
10 10 1 
10 2 10 
10 10 1 
10 
■010100- 



Na 



■0 1 0- 
10 10 
10 10 1 
10 10 
10 

■0 1 0- 



■0 1 0- 
10 
10 
10 
10 

-0 1- 



1' 
10 
10 10 
10 

1 

1 1 0- 
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2 


3 


4 


5 


6 


1 


1 


2 


3 


4 


5 


6 


2 


2 


1+3 


2+4+6 


3+5 


4 


3 


3 


3 


2+4 + 6 


1+3+3+5 


2+4+6 


3 


2+4 


4 


4 


3+5 


2+4+6 


1+3 


2 


3 


5 


5 


4 


3 


2 


1 


6 


6 


6 


3 


2+4 


3 


6 


1 + 5 



Zi = 


XiX*i 


+ X 5 X*5 


+ X 6 X*6, 




Z-i®\ — 


XixX 


+XiX*5 


+x 3 xl 


— ^1(812) 


= 


XzXX 


+x 3 xl 


+ {X2 + Xa)xI 


— ^103' 




XiX*i 


+X2X5 


+X3X6 


— ^1045 


z 5 = 


XsXl 


+XiXt 


+x 6 xt 




z 6 = 


X 6 X*i 


+M*5 


+(xi + x 5 )x* 6 - 





Xi = Xi + Xr, X2 = X2 + X6 + Xs, Xs = 
Xa = X4 + Xe + X10, X5 = X5 + Xn, 



X3 + X5 + X7 + X9, 

Xe = XA + Xs- 



Table 5: The graph fusion algebra of Eq and its twisted partition functions in terms of 
extended characters. The ambichiral subalgebra is shown bold and the extended characters 
are given in terms of the A u characters. 
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For b G G, the rectangular matrix V b = {ni a b )i & A g ^ l ,a€G is called an intertwiner because it 
intertwines the fused adjacency matrices: 

V b = V b m, for all i G A L . (2.12) 

One view of these graph algebras is that nodes label bimodules and edges are homo- 
morphisms between these bimodules. In the case of the graph G, the edges describe the 
homomorphisms arising from tensoring with the fundamental bimodule, the result is isomor- 
phic to the direct sum of the bimodules which are adjacent to it on the graph. For Type I 
models, one can associate a graph G a to each vertex in the same manner, by placing N a b c 
edges between the vertex b and the vertex c. For Type II models, this construction fails. 

The graph fusion algebra of a solvable A-D-E lattice model built on the graph G is 
not the graph fusion algebra itself. Rather, this latter graph encodes the fusion algebra of a 
Wess-Zumino-Witten s£(2) g __ 2 (WZW) theory. The solvable A-D-E lattice model is actually 
associated with a minimal model whose fusion algebra is given by the tensor product graph 
Ag-2 x G where g is the Coxeter number of G. A vertex of this tensor product graph is 
of the form (r, a) G A g _ 2 x G and is adjacent with the vertex (r', b) whenever r and r' are 
neighbours in A g _ 2 and a and b are neighbours in G. When G = A g _i, it is customary 
to denote such a vertex with the Kac labels (r,s). Moreover, this tensor product graph is 
quotiented by the Kac table symmetry 

(r,s) ~ (jg-l- r ,g-s) if G = Vi. (2-13) 
(r, a) ~ (g — 1 — r, cr(a)) otherwise. (2-14) 

where a is the identity for D 2 i, £7, E 8 and Z 2 graph automorphism for D 2 e and E$. The 
graph fusion algebra for this quotient is A g _ 2 <g> Gj ~ and hence it is generated by the two 
s£(2) WZW subalgebras 1 <g> G and A g _ 2 <g> 1. 




Figure 2: The graph {T 2: D^) coding the fusion algebra of the three-state Potts model. The 
solid lines encode the adjacencies of the nodes on the tensor product graph. The nodes 
(r,a) = {(1,1), (1,2), (1,3), (1,4), (2,1), (2, 2), (2, 3), (2, 4)} in (T 2 , D A ) are labelled by the 
associated fields {1, rj, e, £, a, aty respectively. 

For example, the three-state Potts model is the minimal model associated with the 
graph £) 4 . Its fusion algebra is the graph algebra (At, D^)/^ = {T 2 ,D^) pictured in Fig. |21 
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It has eight vertices and two types of edges, corresponding to the action of the generators 
e, 77, of each subalgebra, T 2 = {/, e} and L> 4 = {I, r], ip, where we use the labelling of 
Fig. |21 You read off that £77 = e + a + cfi and £e = £ + 77. Assuming associativity, you can 
then infer for example that = / + e + ?/> + ?/^ + a~ + o^ and work out the rest of the algebra 
table, which is also given by the tensor product of the N matrices for (listed in Table |2J) 
with the fusion matrices of T 2 , namely 

M-P. (2.15) 

2.3 Ocneanu graph algebras and twisted partition functions 

In this section we review the properties of Ocneanu (quantum double) graphs and their fusion 
algebras 0,0,0, 15, 3|. The related Double Triangle Algebra (DTA) governs many aspects 



of both the statistical mechanics models and their associated conformal field theories. It is 
called the algebra of quantum symmetries of the problem. We are interested here in the 
algebra of twisted boundary conditions on the torus, called the twisted fusion algebra 0,0]. 
We discuss the twisted boundary conditions and associated twisted partition functions briefly 
summarizing the work of Petkova and Zuber 0, |^ . We then present an alternative approach 
using tensor products 14, 15, 3|. 



The twisted boundary conditions are encoded by the Ocneanu [14J graph fusion algebra 
G whose structure constants are denoted by N xy z for x,y,z G G. The matrices N y = {N xy z } 
satisfy 

N x N y = J2ti xy z N z . (2.16) 

zed 

Petkova and Zuber give explicit expressions for these structure constants for the A-D-E 
graphs. The seam index 

x= na,b, K )e(H,H,Z2), WZW 
1 (r, a, b, k) G (A g _ 2 , H, H, Z 2 ), minimal 

labels conformal twisted boundary conditions or seams. The index k = 1,2 labels the 
automorphisms 77 = <j K_1 = /, a. The seams x are not all distinct due to quantum symmetry 



x = (r, a, 6, k) = x' = (r, a', 6', k) <^> Z(r,a,6,«) = Z( r t >af>v>K >) (2.18) 

that is, seams giving rise to the same twisted partition functions are considered equivalent. 
In some cases, such as -D 4 , it is necessary to use unspecialized characters to see the full 
quantum symmetry. For the WZW models, it suffices to take x of the form 



S G Ag-!, 


G = 


A 9-l 


(a,K) G (D 2i ,Z 2 ), 


G = 




S G Au-U 


G = 


D 2 e+i 


(a,b) G (G, G), 


G = 


^6,8 


k (a, b) G (D W ,D W ), 


G = 


E 7 



and similarly for the minimal models with r G A g _ 2 added. 
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The modular invariant torus partition functions of the A-D-E minimal models are 



aeT 



-£(<?)= %s),(r',s')Xr,s(9)Xr',s'(9), #(r,«),(r',s') = $r 

(r,s),(r',s') 

where n s = n s (H) are the fusion matrices of the Type I parent H of G 



(2.20) 



T = Ti 



{1,2,. ..,L}, 




G 




{1,3,5,. ..,2* - 


1,2*}, 


G 


= A« 


{1,2,3,. ..,4£- 


1}, 


G 


= £^+1 


{1,5,6}, 




G 


= ^6 


{1,3,5,7,9,10}, 




G 


= £ 7 


{1,7}, 




G 


= E$ 



(2.21) 



and the involutive twist ( is the identity for Type I theories but for Type II theories has the 
action 

s ^ U - s, s = 2, 4, . . . , 2£ - 2, G = D n+ i 
{1, 3, 5, 7, 9, 10} ^ {1, 9, 5, 7, 3, 10}, G = E 7 



(2.22) 



The twisted partition functions are given by the toric matrices 

*W»)(?) = £ N^r 2)r " xw,y(?) x^fo), = E 

(r',s'),(r",s") 



where 



T 2 



{2,4,. ..,2^-2}, G = D 2 , 
Ti, otherwise 



(2.23) 
(2.24) 



except for the special seams denoted x = (r,X) oi E 7 which are given by Z( Ti x)(q) = 

Z( r ,6,2,l)(q) — Z(r,4,2,l)(q)- 

The Ocneanu graph fusion algebra is a double graph algebra combining left and right 
copies of the graph fusion algebra H connected through a left-right symmetric subalgebra 
called the ambichiral subalgebra. For graphs G of Type I, H is G itself. For graphs G of 
Type II, it involves the parent graph algebra H and a twist (. Although the graph algebras 
of Type II theories are not proper graph algebras, the parent graph H is always of Type I. 
The Ocneanu graph has two types of edges (plain and dashed), corresponding to the action 
of the left and right copies of the generator of H. 

An alternative construction of the Ocnean u gr aph fusion algebra, which emphasizes the 



chiral and ambichiral structure, can be given 1J, [l5|, [16J in terms of tensor products over 
the subalgebra T. We start with an example and consider G = -D 4 . The nodes T = {1, 3, 4} 
form a Z 3 subalgebra of the graph fusion algebra G corresponding to the extended chiral 
algebra. This gives rise to the ambichiral subalgebra of G. One can then construct the tensor 
product algebra over T 

D 4 ® z , D 4 (2.25) 

The tensor product <S>t over T means that a (g> tb = at <g> b for all a, b e G and t G T. In 
D4 for example 2cg>3 = 23(g>l = 2cg>l. So this algebra has six distinct elements, 1<8>1, 



11 




Figure 3: The Ocneanu graph D A . The nodes can be labelled by x = (a, b) or (a, k). We have 
set a = a<g>l, 2' = 1(g) 2 with 2®2 decomposed into l',3',4'. The nodes {1, 2, 3, 4, 1', 2', 3', 4'} 
are alternatively labelled by (a, k) = {(1, 1), (2, 1), (3, 1), (4, 1), (1, 2), (2, 2), (3, 2), (4, 2)} 
respectively. 




6' 



Figure 4: The Ocneanu graph D 2 e for i = 3. 
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1®2 = 3®2 = 4®2, 1®3 = 3® 1, 1 ® 4 = 4 ® 1 ,2®1 = 2®3 = 2®4, and 2 ® 2. It 
clearly has two generators, a left and a right, Hl — 2® 1 and D R = 1 ® 2. They generate the 
ie/t respectively ng/it chiral subalgebras. One would like to encode this algebra in a graph 
as previously but there is an obstruction: 

(2 ® 2)(1 ® 2) = 2 ® (1 + 3 + 4) = 2 ® 1 + 2 ® 1 + 2 ® 1 (2.26) 

so while there is only one edge from 2 ® 1 to 2 ® 2 = (2 ® 1)(1 ® 2), there would be three in 
the opposite direction. This problem is solved by splitting the node 2 ® 2 into three different 
ones, 1', 3', 4', using a non central extension of the algebra by an algebra of 2 x 2-matrices [lfij ]. 
The detail of this extension is not needed to compute twisted partition functions but we note 
that it leads to non-commutativity. One then obtains the graph D 4 presented in Fig. El 

The same procedure works for G — H — D 2 g (see Fig. HJ). It has a Je + i subalgebra 
generated by the odd vertices (the forked extremities are both taken as odd) over which the 
tensor square is taken and extended by an algebra of 2 x 2-matrices. The result is encoded 
in a graph D 2i , depicted in Fig. HI where the ambichiral 1, 3, . . . , 2£ — 3, 2£ — 1, 2£ nodes are 
duplicated by 1', 3', . . . , (2£ — 3)', (2£ — 1)', (2£)', and where the even vertices have a left and 
aright counterpart 2® 1, 4® 1, . . . , 2(£- 1) ® 1 and 1 ® 2, 1 ® 4, . . . , 1 ® 2(1- 1). The algebra 
structure can be worked out beginning with 1', which satisfies 

l'(l®2) = 2®1, (2.27) 
1'(2®1) = 1®2, 

while 2 ® 2 = V + 3'. 

The Type I exceptional cases E§ and E$ are simpler as there is no need to extend the 
algebra. The subalgebras in these cases are generated by T = {1, 5, 6} and {1, 7} respectively. 
The Al case is even simpler as we tensorise over the full graph algebra, yielding the same 
algebra back again Al ®a l Al ~ Al so all the elements are ambichiral. 

For the Type II models, G = D 2 e+i,E 7 , the Ocneanu graph algebra is defined through 
the square tensor of the parent theory H = Au-i,Dxq twisted by the involution (. For 
G = D 2i+1 , the twisted fusion algebra is defined by A^i ®^ A^_i where a ® b = a((b) ® 1 = 
1 ® ((a)b. It is the graph algebra of the Ocneanu graph D 2 £ + i defined by 4£— 1 vertices 
forming two A^_i graphs with plain and dashed edges, sharing the same numbering for 
the odd vertices (forming Exp(.D 2 £+i)) but where the even ones are flipped, as pictured in 
Fig. El For G = Ef, the parent theory is H = Diq and the automorphism is given by 
interchanging the nodes 3 and 9, so that the twist fusion algebra is D w ®t,( D w where T is 
the .Dio ambichiral subalgebra (its odd vertices, counting forked vertices as both odd). The 
left (right) chiral subalgebra is generated by the left (right) generator and the ambichiral 
subalgebra. They are both isomorphic to the primitive graph algebra H and, as in the Type 
I case, the ambichiral subalgebra is their intersection. 

To summarize, the vertices of the Ocneanu graphs are given by distinct pairs (a, b) G H 
(H is the graph G or its parent graph), coding a left and a right element, or equivalently 
for the D 2 i models, a pair (a,n) with k G {1,2}. The twist fusion algebras just described 
are for the s£(2) 9 _ 2 WZW theories. The Ocneanu graph fusion algebras associated with 
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1 




Figure 5: The Ocneanu graph D 2 e+i- 

A-D-E minimal models are more involved. These are given by two copies of the graph 
algebra (A g _ 2 , H) / ~ whose vertices are pairs (r,a) with r G Ah_ 2 and a G H, subject to 
the Kac table symmetry (r, s) ~ (g — 1 — r, g — s) if H = A g ^\ and (r, a) ~ {g — 1 — r, a) 
otherwise. The vertices are, in general, labelled by x = (r,a,b,K,) G {A g _2, H, H, Z 2 ). In 
the (A 4 , D4) case for example, the Ocneanu graph is simply the tensor product of D4, the 
Ocneanu graph of D4, with the tadpole T 2 . It has three types of edges, associated with the 
action of the left, the right and the tadpole generators. 

2.4 Twisted partition functions and extended characters 

Given an A-D-E graph G, we associate to each node a G G, an extended character x a {q) 
which is a generating function in the modular parameter q. The ordinary characters Xi{q) 
are the A-type characters. The extended characters are 

Xa {q)= Yl n n a Xi(q) (2-28) 

Explicitly, in the L> 4 case, 

Xl = Xl+ X5, X2 = X2 + X4, X3 = X4 = X3- ( 2 -29) 

We use the extended characters mostly for Type I graphs, that is, either G itself or its parent 
H if G is of Type II. 

The twisted partition functions, which are sesquilinear combination of characters, can 
now be written in terms of extended characters in the following way jHSBEB]- First, the 
modular invariant partition function corresponding to the unit vertex of the Ocneanu graph 
is associated with a sum over the subalgebra T 

Z{q):=Y,Uq){xaa)W ■ (2-30) 
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where the automorphism ( is the identity for the Type I models. For Type II theories, the 
sum over T is a sum over (a subset of) the parent graph. In the D 2 £+i case for example, 
we sum over Au-i (they are all ambichiral). Moreover, the vertex a, associated with Xa{c[), 
is paired by the twist ( with Xc( a ) (<?)*■ It is also a sesquilinear form in terms of ordinary 
characters and its diagonal part gives the exponents Exp(G) of Table ^ This partition 
function has the property of modular invariance in the modular parameter q = exp(27rir), 
that is to say, Z(l + r) = Zir) and Z(—1/t) = Zir). The other twisted partition functions 
are not modular invariant. 

The twisted partition function of other elements a ®t b of the Ocneanu graph, with or 
without an automorphism seam rj = a K ~ x , are obtained by the action of this element on the 
terms of Z(q) 

zitM ■= E (xc(c)*(g))* = E N ac d N Hi cfUQ) (Uq)Y (2.3i) 

c£T K c£T K ,d,eeH 

where ac and cb denote the product in the graph fusion algebra, that is to say, it is a 
linear combination of elements of the algebra and the characters are extended by linearity 
in their indices. The graph H = G for Type I theories but is the parent graph for Type 
II theories. Petkova and Zuber p indicated that the twisted partition functions could be 
written in terms of the extended characters but did not write down an explicit formula. 
Notice that this implies that the N structure constants involved in the calculation are all 
positive integers since no use is made of the graph fusion algebra for Type II models, only 
for the parent model which is always of Type I. For example, the action of the left generator 
changes the left term in x a (q) (xc(a)(?)) to a sum over the neighbours of a in H 

Zn L (q) = Z 2 ®i(q) = E X2a(q) (x«a)(?))* = E Xb(q) (&(«)(?))* ■ 

Notice that Z a( ^ T f, = {Z^ Ta )* and, in particular for the Type I theories, the ambichiral 
twisted partition functions are real. A complete list of these twisted partition functions in 
terms of ordinary characters is given in N, ha. For the cases D e and E 6 we give a list of 
them in terms of extended characters in Tables 0] and El 

The twisted partition functions just described are for the si(2) g _ 2 , or WZW models. The 
twisted partition functions for the minimal models (J2.23|) involve pairs (r, a) G A g _ 2 x H/ ~ 
of indices in place of single indices a G H . Here, g is the Coxeter number of H and ~ is 
the Kac table symmetry. The ambichiral subalgebra is the product of the two ambichiral 
subalgebras for each graph. In the case (A 4 ,_D 4 )/ ~= (T 2 ,-D 4 ), for example, the modular 
invariant partition function is 

Zi = E XxXl = E = Xl,i + Xl,i + Xls + Xl, 3 + Xi, 4 + X3.4 (2-32) 

x£T re{l,3},ae{l,3,4}cD 4 

where we have kept the labels 1 = 4 and 2 = 3 for r G T 2 = A4/Z2. The other 15 twisted 
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partition functions are obtained by action of the twisted fusion algebra. For example, 

Z (l,2)tg>l = Xl,2(Xl,l + Xl,3 + Xl,4)*+X3 ! 2(X3,l+X3,3+X3,4)* (2.33) 
2(3,1)91 = X3,lX*,l + (Xl,l + X3,l)X3,l + X3,3Xl,3 + (Xl,3 + X3,3)X3,3 

+X3AX*1,4 + (Xl,4 + X3,4)X3,4 ( 2 - 34 ) 

and so on as listed in Table El 



Zl = Xl,lXl,l + X 3 ,lX3,l + Xl,3Xl,3 + X 3 ,3X3,3 + Xl,4Xl,4 + £3,4X3,4 

z m i = xi,2(xi,i + Xi, 3 + Xi, 4 )* + X3,2(x 3 ,i + x 3 ,3 + X3 A T = z\ m 

Z i> = Xl, 3 Xl,l + £3,3X3,1 + Xl,4Xl,3 + X 3 ,4X3,3 + £l,l£l,4 + £3,1X3,4 = Z ^ 
Z e = £ 3 ,lXl,l + (£l,l + X 3 ,l)X3,l + X 3 ,3Xl,3 + (£s,3 + Xl, 3 )X3,3 + 
X 3 ,4Xl,4 + (£ 3 ,4 + Xl, 4 )X3,4 
Z m = X3 ) 2(Xl4 + Xl,3 + Xl,4)* + (Xl,2 + X3,2)(X3,l+X3 1 3+X3,4)* = Z * m 
Za = £l, 3 £l,l + £3,3X3,1 + Xl,4Xl,3 + £3,4X3,3 + £l,l£l,4 + £ 3 ,lX3,4 

= z a1 

Z\> = \Xi,2\ 2 + 1X3,2 1 2 = Z^> = 

Z e' = Xl,2X3,2 + £3,2X1,2 + IX 3 , 2 | 2 = Z a > = Z aV 

Xr,l ■= Xr,l + Xr,5, Xr,2 ■= Xr,2 + Xr,4, £r,3 = £r,4 ■= Xr,3- 

Table 6: The (At, D4) twisted partition functions of the 3-state Potts model. The extended 
characters are given in terms of the Virasoro minimal characters Xr,s- 



3 Lattice Realization of Twisted Boundary Conditions 



3.1 A-D-E lattice models 

A solvable 0] A-D-E lattice model 0, 18] is associated with a graph G, of A, D or E 
type. We place spins on the sites of the square lattice, where the spin states are taken to 
be the nodes of the graph G and neighbouring sites on the lattice must be neighbouring 
nodes on the graph. The probability distribution of spins is defined by the critical (unfused) 
Boltzmann weight of each face (or plaquette) of spins, depending on a spectral parameter u. 
For four spins a,b,c,d G G such that (a, b), (b, c), (c, d), (d, a) are pairs of neighbours in G, 
the Boltzmann weight is 



W 



11 



u 



11 



s(X - u) 5 ac + s(u) 



1pb1pd 



>bd 



(3.1) 



and zero otherwise. Here, g is the Coxeter number of G, A = ~, s(u) = and ip a is 

the entry, associated with the node a, of the Perron-Frobenius eigenvector of the adjacency 
matrix G. 
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These Boltzmann weights can be represented Q by a local face operator Xj\ 



u 



Xj(u) 




s(X — u)l + s(u)ej 



(3.2) 



in the Temperley-Lieb algebra T(N, X) 

s(2A) e,- 



6 j 6fc 6 j 
6 j 6fc 



6fc Cj, 



\j -k\=l 
\j -k\>l 



(3.3) 



where ej = Xj(X) is a Temperley-Lieb generator and j is an integer labelling the iV positions 
along a row of the lattice. 



3.2 Fusion projectors 

Each of the A-D-E models gives rise to a hierarchy of fused models whose Boltzmann weights 
we are going to describe. They are associated with blocks of faces where the internal spins 
are summed over in a particular way. 

We first define recursively the fusion operators PJ, for r — 1,2, . . . , g as follows 



^ -PJ+i 1 Xj(—(r-2)X) P T j~l , r>3 



(3.4) 



where = s(k\) and j is appropriately restricted |4|. Thus, PJ can be expressed as a 
function of ej, ej + i, . . . , e 3 - + ( r _3). In particular, 




S 2 



(3.5) 



We shall represent the fusion operators diagrammatically as 



pr 

j 



j'-i j' 



(3.6) 



j+r— '3 j'+r— 2 



It is easy to show that this operator is in fact a projector. Moreover, 



pr pr 



pr pr' _ pr 
^3 ~ ^' ' 



for < j' — j < r — r' . 



(3.7) 
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In particular, the local face operator = ls a projector orthogonal to all the PJ 

for < j — j' < r — 3. This fact is a defining property of the orthogonal projector PJ 

j+r-3 

ImPJ= p| Kere fc . (3.8) 

k=j 

Clearly, we can decompose the projector PJ onto the space of paths with given end points: 
P r (a,b) is the fusion projector acting on paths from a to b in r — 1 steps. Its rank is given 
by the fused adjacency matrix entries 

Rank (P r (a, 6) ) = n ra b . (3.9) 

The +1 eigenvectors of P r (a, b) are thus indexed by an integer 7 = 1, 2, . . . , n r a b referred 
to as the bond variable. We denote these eigenvectors by U^a, b) and call them fusion 
vectors or essential paths. Explicitly, in the representation (|3.1|) of the Temperley-Lieb 
algebra T(r — 2, A), these generators act on the paths from a to b in r — 1 steps as 

e/c(ao, ai, . . . , a^-i, a^, a^+i, • • • , a r -2, b) = (3.10) 
o a ^_ liahH ~Jv2 — 7^5— l a o, a i, • • • , c, a fc+ i, . . . , a r _ 2 , oj. 

In the D and E cases multiplicities occur and there is some freedom in the choice of fusion 
vectors corresponding to a unitary change of basis. In the Ai case, however, there is a unique 
fusion path. As an example, there are two paths on A^, going from the node 2 to itself in 
2 steps, namely (2,1,2) and (2,3,2). As they both backtrack, the fusion vector C/^(2,2) 
is unique, proportional to their difference '0^ 2 (2, 1, 2) — V ; i /2 (2, 3, 2) and the fused adjacency 
matrix entry is n 32 2 = 1. In the D4 case, there is the path (2,4,2) as well, so that there 
are two linearly independent fusion vectors, the previous one and (2,3,2) — (2,4,2) or any 
similar linear combination (notice that ip3 = ^4). The fused adjacency matrix has a two as 
the corresponding entry n 32 2 = 2. The general form of the unique A^ fusion vector at fusion 
level s < L between the vertex 2 and the vertex s — 1 is given by the following formula of 
cancelling alternating backtracking paths, generalizing the one just described for s = 3 

175(2, s - 1) = (V>2 1/2 (2, 1, 2, 3, . . . , s-1) - 3, 2, 3, . . . , s-l) 

+ ^7(2,3,4,3,..., S -l) + ... + (-l) s ^; 1/2 (2,...,s-l, S ,s-l)) (3.11) 

and similarly for Dg with fusion level s < i — 1. But at fusion level s = I — 1, the fork gives 
rise to a two dimensional space of fusion vectors. One choice of orthonormal basis is given 
(with v the appropriate normalisation constant) by 



c/ £ r 1 (2, t-i) = \ {$T$T% 1, 2, 3, . . . , 1-2) - ^s 1/2 (2, 3, 2, 3, 



'-2) 



+ Vr>4 (2,3, 4, 3,..., £-2) + ... (3.12) 

+ 5(-i)*-V3Vi?(2, 2, i- 1, t - 2) + ^(-1)^7^(2, • • • , e-2, e, £-2] 

U e f\2, £-2) = ^= ((2, ...,£-2,£-l,£-2)-(2,...,£-2,£,£-2) 
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3.3 Fused face operators 

The fusion projectors allow us to define the (p, g)-fused face operators consisting of q rows 
of p local face operators with relative shifts in the spectral parameter by ±A from one face 
to the next 

u+(q-l)\ 
«+(<MP)A 



j+p-1 



Xf{u) 




j-1 j+q-1 j+p+q-2 




(3.13) 



u— (p-l)A 

The position of the projectors and spectral parameters can be altered by pushing-through 



u-(p-l)X 






u+(q-p)X u+(q-l)X 

These properties imply several others, namely the Transposition Symmetry 

X™(u) T = Xf(u+(q-p)\) 
the Generalized Yang-Baxter Equation (GYBE) 



(3.14) 



(3.15) 




the Inversion Relation 



X™{u)X?{-u) 





(3.16) 



p-1 q~l 

where s P9 (u) = J j j^J s(u + (i—j+k)\) (we will also use the notation for q — 1) and the 

j=0 fc=0 

Abelian Property 

Xf{u + {p-l)\)Xf{v + fo-l) A) = *f (t; + (p — l)X)X QP (u + (g-l)A) . 
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The braid limits of the fused face operators are 



- (g+l)P9 - 



Xf(ieoo) 



lim 

«— >ieoo 



X pq (u) 



±1 



It follows that 



X? (+zoo) = Xf (-ioo)* 
and the inversion relation becomes 

Xf (zeoo)Xf (ieoo)t = P/ +1 Pjj; 



(3.18) 



(3.19) 



(3.20) 



so that these operators are unitary. 

The fused face operators, contracted against the fusion vectors, yield the (p, g)-fused 
Boltzmann weights which depend not only on the spins at the four corners but also on the 
bond variables on the edges 

d U 1 (d,c) c 

d 7 c 



W pg \ S (3 
a a b 



u 



u /p 



u 




&TW (3.21) 



a TT p ~ 1 r 



where the function Sq' 3 ^ 1 (m) eliminates some scalar factors common to all the spin configura- 
tions which appear in the process of fusion. In the Al case the bond variables are redundant. 
The fused Boltzmann weights satisfy the reflection symmetry 



/ d 


7 


c 








(3 


■) 


r 

y a 


a 


b 





Sg^p (it) 



d 5 a 



u 



W qp \ 7 



a 



c 13 b 



u + (q-p)X 



(3.22) 



and Crossing Symmetry 
W pq 




' j'ai'c S^jX-u) 



u 





( a 


s 


d 




w qp \ 


a 




7 


A-uj 




Kb 




c 





(3.23) 



We use these fused Boltzmann weights to construct commuting transfer matrices with seams. 



3.4 Integrable seams on the torus 

Simple integrable seams are modified faces. Surprisingly, they produce some new twisted 
boundary conditions even for the Ising model 19, |2(|. They come in four different types, 
r, s, a and //-type where r G A g _ 2 , s G a E H and G T. Here H = G for Type I 

theories and is the parent for Type II theories. A composite integrable seam x consists of 
several simple seams glued together with four spins c, d, e, / at the corners. The integrable 
seams x = (r, a, b, <j K ~ l ) give rise to the conformal seams (r, a, b, k) in the continuum scaling 
limit. 
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First, let us consider the case r = 1 equivalent to the WZW factor in the minimal coset. 
The adjacency matrix of an integrable seam x = (a, b, <r K ~ 1 ) is given by 



N a N b a k ~\ G of Type I 

n x = { nf H) nf H) a k -\ G of Type II and x ^ X (3.24) 

nf H \nf H) -nf H) ), G = E 7 and x = X 



where n^ H ^ with a G H are the intertwiners [3J of G relative to H 

\\H)m 



T )c = £ %r4 G)m 4 G)m * (3.25) 



meExp(G) ™1 



and X is the special node [8( of the E? Ocneanu graph. The matrices (n x ) c d and (n x )f e 
encode the allowed adjacency between spins c, d at the bottom and /, e at the top of a 
composite seam. Although our interpretation of these matrices is different, these matrices 
coincide exactly with the n x matrices of Petkova and Zuber. Our definition, however, is 
intrinsic to the seam so we do not need to invoke boundary conditions on the cylinder. The 
matrices n x form a (non-faithful) representation of the Ocneanu graph fusion algebra 



Tl X Tly 



Y,N xy z n z (3.26) 



where the explicit expressions for iV are given in 

For the unitary minimal models (A g -2, G), the Ocneanu graphs are the quotient of the 
tensor product graph, A,_ 2 <g> Gf ~, where ~ is given by the Kac table symmetry (j2.14|) . 
There are (g — 2)\G\/2 distinct nodes on this Ocneanu graph and the corresponding graph 
structure constants Nh, } y) = {^(p',y')(p,y)^ p "' y "' > } are given by 

AT {P",y") — n (r,x) n (r',x') n (r",x") , r (A,- 2 )r" ( <j \ x" (1 97\ 

where (p,y) = (p,a,b, re), 

Q(p,y)^ ' ^ 3p,r$y,x ~\~ ^p,g— 1— r ^y,T}(x) (3.28) 

where the indices (p, y) range over 



1<P<^, l<y<\G\, seven, 
l<P<2-2, l<y< 9 -^, G = A g - 1 ,g odd, 



(3.29) 



and the matrices Q act on each index of the tensor product Nr Ag 2 ^ <g> N x to quotient out 
the Kac table symmetry. The matrices 

n M = Q { pJ r ' x) ® N x (3.30) 

now form a (non-faithful) representation of the Ocneanu graph fusion algebra for minimal 
models. 



21 



We construct composite seams of type x = (r,s,rj) and first define W? rll s, the r-type 
seam for the (p, g)-fused model. It is a usual (r — l,g)-fused face (it doesn't depend on 
the horizontal fusion level p) with an extra parameter £ acting as a shift in the spectral 
parameter, and another choice for the removal of common scalar factors 
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_r— 1 q—1 
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An s-type seam with s G is the normalized braid limit of an r-type seam, it does not 
depend on the spectral parameter 




d j c 



(l,s)/i 



:(9+l)(»-l)g- 



lim 



Id 7 
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\ a a 


b 





(3.32) 



In general the s-type weights are complex. The complex conjugate gives the weights in the 
other braid limit £ — > —ioo. By the reflection and crossing symmetries ()3.22|) and ()3.23|) 



d 7 c 







a a b 



# *3 



(l,*,l) 




(3.33) 



These braid-limit face weights provide us with a representation of the braid group. 



It is known 
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22| that discrete symmetries play an important role in twisted boundary 



conditions. In fact, there is an integrable seam corresponding to each discrete symmetry. 
Specifically, the graph automorphisms r\ G T, satisfying G a ^ = G^a)^)-, leave the face 
weights invariant 
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\ a 
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b 





u /i 



w pq 



and act through the special seam [2 



u /) 



77(a) 



r,(6) 




ObT)(a)Vcr)(d) 



u r) ft 



Tj(d) 

! 

r?(a) 



/ 77(d) 7 77(c) 
y 77(a) a 77(6) 



1, %ia d 7^ 0, a = /3 
6 = 77(a), c = 77(d), 
0, otherwise. 



(3.34) 



(3.35) 



Notice that the (r, s,t/) = (1, 1, 1) seam, where 77 = 1 denotes the identity automorphism, is 
the empty seam 

d c 

5 a b 5 cd 5 a/3 n q+lb c . (3.36) 



a (3 
a b 



The push-through property is also satisfied for an 77-type seam. 
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3.5 Construction of a-type seams for Type I theories 

In this section we construct a new type of fusion for Type I models labelled by the nodes 
aeG. The fused Boltzmann weights giving the (1, a)-seams are obtained in the braid 
limits u — > ieoo with e = ±1 and are independent of the spectral parameter u. This new 
type of fusion is associated to the N a graph fusion matrices in exactly the same way as 
the usual fusions are associated to the fused adjacency matrices n s . For Type II models, 
the construction is applied to the parent Type I graphs in the next section. Previously, all 
known fusions were labelled by Young tableaux. In our case, r and s label s£(2) tableaux 
with r or s boxes in a single row. The construction of fusions labelled by nodes of graphs 
is an important step in understanding the graph fusion algebras associated with integrable 
and conformal seams. 

In constructing the a-type seams in this section we focus on vertically unfused face 
weights with q — 1. Recall that the s-type seams are the braid limit of the (s — 1) fused 
Boltzmann weights which are obtained by acting with the fusion projector P? on the face 
operators Xj(u). From (j3.9|) we see that the admissible spins at the corners of the (s — 1) 
fused face weights and the number of horizontal bond variables are given by the non-zero 
entries in the fused adjacency matrices n s with s G ^4 9 -i- From Section l2~2| these fused 
adjacency matrices are linear combinations of the graph fusion matrices 



$>.i a iV a (3.37) 



n s = 

aeG 

This motivates us to define a new type of fusion projector Pj S ^ a (b, c) associated with N a and 
acting on paths between b and c by orthonormally decomposing the s-type projector Pj{b, c) 

P;{b,c)=Y j n s i a Pf a {b,c) (3.38) 

aeG 

where each of the N a fusion projectors Pj s ^ a (6, c) separately satisfies the push-through prop- 
erty (J3.14j) in the braid limit. These projectors P^ a (b,c) which may have complex entries 
are required to satisfy 

1. Decomposition: P s (b, c) = ^2n sl a P {s)a (b, c) (3.39) 

aeG 

2. Orthogonality: P {s) a (6, c) P (s) f (b,c) = 0, a^f, (3.40) 

3. Projection: P (s)a (6, c) 2 = P {s)a {b, c), P (s)a {b, c) f = P (s)a (6, c) (3.41) 

4. Adjacency: Rank(P (s)a (6, c)) = N ab c , n sl a ^ 0. (3.42) 

By the above conditions (|3.39|) - (|3.42|) . P s (b,c) and {P^ a (b, c)} a( z G are simultaneously 
diagonalizable with a common set of eigenvectors. It follows that P^^o, c) can be decom- 
posed as 

P^(b,c) = ±U^(b,c)U^ a (b,cy (3.43) 
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where U ' ^ (5, c) are the simultaneous eigenvectors of P s (b, c) and {P^ a (b, c)} a6 c satisfying 



P°(b, c) uf\b, c) = PW»(6, c) U { ^\b, c) = U ( ; )a (b, c) 
P^f(b,c)U ( ;' )a (b,c)=0, a^f. 



(3.44) 



Thus, the construction of the new type of fusion at u — > ieoo, with e = ±1, is equivalent 
to finding the appropriate ortho normal basis for the fusion vectors {U^l {b, c)} of P s (b,c) 
satisfying the conditions 



r yWf t 



a [s) a == 

U efl (b,c) 



for all fi,u,a^ f, N ab c , N fe d ^ 



which follows from the push-through property and the orthogonality of P^ a ()3.4())) . Note 
that the N a fusion vectors at u — > ieoo are different for e = +1 and e = — 1. By (j3.19|) . 
t// (6, c)* is the iV a fusion vector at u — > —ieoo. Hence, the iV a fusion projectors at 



(3.45) 



K 



±ioo are complex conjugates 



P {s)a . (b c) = P 



[s)a 
u— >+ioo 



(6,c)*. 



(3.46) 



From now on, unless otherwise stated, the terms N a fusion projector and N a fusion vector 
refer to the braid limit u — > -Moo. We emphasize that p^ a does not satisfy the push- 
through property if the face operator X^~ 1 ^ 1 (m) depends on u, that is, if we move away from 
the braid limit. 

The N a fusion vectors can be obtained by solving ([3.45)1 . However, there is a more 
convenient approach taking advantage of the unitarity of the face operators in the braid 
limit (|3.20jl . This provides a unitary transformation between essential paths, so that the 
unknown N a fusion vectors can be obtained from known ones by a unitary transformation. 

The fusion vectors U s (e, d) with N ae d = 1 are automatically the same as the N a fusion 
vectors due to the adjacency condition ()3.42|) 



U s (e,d) = U is)a (e,d), N a 



(3.47) 



r( S )°/ 



So suppose that U (e, d) is known for some given s, d, e. Then for any b,c G G with 
N ab c ^ satisfying 

n 2 d 7^ and n 2e b ' n sb > c = n sb c 8 Vjb (3.48) 



(s)a 



the N a fusion vector U (6, c) is given by the unitary transformation 



U„ (e, d) i * J7 (6, c) : J7 (6, c) 




(3.49) 
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Clearly, it is a fusion vector of P s (b, c) by the push-through property of the ordinary projec- 
tor. It is also an N a fusion vector because 



e U„ (e,d') d' 



Tyw a 
U„ (6,c) 




C7„ (e, d)T . J7 ( 



e, d) 5, 



o a ,f <W 6d' } d 



(3.50) 

by the inversion relation ()3.2()j) . which is guaranteed to be valid by ()3.48j) since, for any V 

adjacent to e, we have Rank(P^^ a (6', c)) = n s b c 5b>,b- Thus, the N a fusion vectors tj\J (b,c) 
are determined uniquely up to a phase factor and unitary gauge transformations within 
blocks given by N ab c > 1. To be consistent, one should check that ()3.49|) satisfies ()3.45|) for 
all the other admissible spins. By the crossing symmetry ()3.23|) and ([3. 19)1 . (|3.49|) can be 
expressed as the linear combination of ordinary fusion vectors 



0?"(M = E 



i>di>b „, (bye 



We now illustrate the calculations of N a fusion vectors for the case G — Eq. In Eq, 
n s — N s for 1 < s < 3, so it remains to find the N a fusion vectors for N^, N 5 and N$. First, 
n 4 can be decomposed into iV 4 and N 6 



n>4 





iV 4 + iV 6 (3.52) 



and only the N a fusion vectors for (6, c) = (2, 3), (3, 2), (3, 4), (4, 3) need to be determined 
that 



* (4) 4 ~ (A) 4 

Consider (6, c) = (2,3). By (EOHJ) . CT ; (1,4) ^ IT ; (2,3) so it follows from (J3~5TJ) 



^ft'HE^sM? i O^ 2 ' 3 ' (3 ' 53) 

and similarly t/ (4)6 (l,6) ^ L/ (4)6 (2,3). For (6, c) = (3,2), {C/ (4)4 (2, 3), L/ (4)6 (2, 3)} i-> 
{IT ; (3,2),ET ; (3,2)}, so CT ' (3,2) is obtained by applying fl3~H to ipH)5|) 



where 



(1,4,1) 



3 7 2 
2 4 3 



) = ^C/ 4 Q (2,3)t. £ / (4)4 (2,3) W (W Y 



3 7 2 
2 a 3 



iy (1 



(1,4,1) 



2 a 3 
1 1 4 



372 



(w) V 2 a 3 



(3.54) 

(3.55) 
(3.56) 
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is the s-type seam with the bond variable 4 specifically chosen so that the corresponding 
fusion vector between (2,3) is N4. Similarly we obtain ^ 6 (3, 2). For the other cases, we 
apply the same procedure 

{t/ W1 (6,3),£/ (4)6 (6,5)} M {t/ W4 (3,4),t/ W6 (3,4)} 
{ ^» 1 (5,2),t/ W6 (5,6)}^ { ^ )4 (4,3),(7 W6 (4,3)} ^ 

Finally, we repeat the whole process for 115 = N3 + N$ to obtain the N$ fusion vectors. 

In summary, the N a fusion vectors for E§ derived from and are listed below in 
terms of the given path basis with A = n/12. 
For n 4 = N A + N 6 : 

(2, 3) Basis: {(2, 1, 2, 3), (2, 3, 2, 3), (2, 3, 4, 3), (2, 3, 6, 3)} 

4 (2j 3) = ^zi f( ^3 + 1) V2 e -f *, - e -^, e ¥*. -(^3 + l)^e^) (3 . 58) 
C/(4)6 (2)3) = v^i ^/2e-^i ) _^_ i)i/2 e -^,_(V3 + i)i/2 e ^,V2e?*) 

(3, 2) Basis: {(3, 2, 1, 2), (3, 2, 3, 2), (3, 4, 3, 2), (3, 6, 3, 2)} 

j7(4)4 (3>2)= v^l (_(^3+ l)l/2 e 5Ai ) e 5A t) _ e -A l)( ^3 + 1 )l/2 e -4Ai) (3.59) 
#(4)6(3,2) = y^(_^ e 5Ai )( ^3_ 1) l/2 e 5Ai )( ^3 +1) l/2 e -4Ai ) _^ 2e -2Ai) 



(3,4) Basis: {(3, 2, 3, 4), (3, 4, 3, 4), (3, 6, 3, 4), (3, 4, 5, 4)} 

[>(4) 4 (3) 4) = 4=1 ^-^ _ e ^ ; + i)i/2 e -^< ) (^3 + i)i/2 e ^ 

C/(4) 6 (3) 4) = y^i ^_ ( ^ + 1 )i/2 e -i^i > _ i)i/2 e ^< > V2e"^ 1 , ^t*) 



(3.60) 



(4,3) Basis: {(4, 3, 2, 3), (4, 3, 4, 3), (4, 5, 4, 3), (4, 3, 6, 3)} 

^(4) 4 (4) 3) = y^i 4j _ e -^i > (^3 + 1) i/2 e -^i j + i)i/2 e ^i) (3 61) 

[>(4)6 (4) 3) = 4_i ^_ ( ^ + ^1/2^^ _ !)i/2 e -f i ? > /2 e -^< ) ^2 e f ^ 

For n 5 = N 3 + N 5 : 

(2, 4) Basis: {(2, 1, 2, 3, 4), (2, 3, 2, 3, 4), (2, 3, 4, 3, 4), (2, 3, 4, 5, 4), (2, 3, 6, 3, 4)} 

#(5) 5( 2j 4) = 4=1 (-y^e 3 ^, (V3 - l) 1 ^ 3 * (^3 _ l)l/2 e -3Ai j _ v ^ e "3Ai ) (3 62 ) 

£/( 5 ) 3 (2, 4) = (-^ - 1)V2 (- ^2e 3A % (V3 - l) 1 /^* _(^3 _ i)V* e - Xi , ^g-Ai, e" 5Ai ) 



(3, 3) Basis: {(3, 2, 1, 2, 3), (3, 2, 3, 2, 3), (3, 2, 3, 4, 3), (3, 2, 3, 6, 3), (3, 4, 3, 2, 3), (3, 4, 3, 4, 3) 
(3, 4, 3, 6, 3), (3, 6, 3, 2, 3), (3, 6, 3, 4, 3), (3, 6, 3, 6, 3), (3, 4, 5, 4, 3)} 

#(5) 5( 3j 3) = (4=1)3/2 ((^3 + 1 )l/2 j _ 1; ^ (^3 + l)l/2 e -3Ai > _^ X _(^ + 1)1/2^3* 

(V3 + l) 1 /2 e 3A< j _(^3 + l)l/2 e 3Ai > Q) _(^ + 1} l/2) 

[/f ) 3 ( 3 , 3) = 1(1 - ^)V2 ((^3 + l) 1 /2 e -2Ai ) _ e -2Ai ; e -2A l) Q) e _2Ai ) _ e -2A i; Qj Q> Q) 

0,(V3 + l) 1/2 e- 2Ai ) 

f/f 3 (3, 3) = 3-1 (%i) 3 / 2 ((V3 + l) 1 /2 e ~5Ai ) _ e -5A 4) _^ 3e Ai ) (2(^3 + 1 ))l/2 e -Ai > 

> /3 e Ai > e -5Ai } _( 2 (^3 + l)) 1 /2 e -Ai ) _(2(V3 + l)) 1 / 2 e 3Ai , 
(2( v / 3 + l)) 1 /2 e 3Ai ) 0, -(V3 + l) 1/2 e" 5Ai ) 



(3.63) 
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or any other orthonormal basis spanned by C/} 5 ^ 3 (3, 3) and L^ ^ °(3, 3). 

(4, 2) Basis: {(4, 3, 2, 1, 2), (4, 3, 2, 3, 2), (4, 3, 4, 3, 2), (4, 5, 4, 3, 2), (4, 3, 6, 3, 2)} 

U® 5( 4j 2 ) = V3pl (_^2e-3Ai ; (^3 _ IJl/Ze-Stt (% /3 _ l)l/2 e 3Ai 5 _^2 e 3Ai ; _^) (3 _g 4) 
#(5) 3( 4j 2 ) = (_| _ 1)1/2 ( V ^ e -Ai j _(^3 _ l)l/2 e -Ai ) (^3 _ l)l/2 e 3W _^2 g 3Ai ? e -5Ai) 

The a-type seam weights, associated with the N a fusion projectors, are given by 



r(5)3. 



^(1,0,1) 





V 


d 




(l,a) 








) = 




/' 


c j 




N, 



j (3 + l)(s~l) \ 

e 2 
lim 

it-noc s(w + £) s!T 1 2 (m + ^) 



~ (s)a t 

17„ M) 



(3.65) 



r "y(s) a 



where s is chosen so that the corresponding fused adjacency matrix contains N a (12.11)1 . Note 
that the choice of s to obtain a particular a-type seam is not unique. In the Eq case, for 
example, we can obtain the Na seam from 714 = N4 + Nq or from Uq = N 2 + Na. Remarkably, 
the weights for such cases agree up to a unitary similarity transformation and give the same 
spectra for the transfer matrices. 

As an example we give the a-type seam weights explicitly for Eq: 



N 



n^ 
n-n 



e 2 , 



-2iA 



-e 2 , 



n-n- 
n^n 



e 2 



iA 



. e -|iA 



n 



_ e -|iA 



9 - -1+V3 p 3i\ 



N R : 



n 



2. 3 



— e2 



\ 3 



3 4 



— e 2 



3 — 4 



n-n 



6 6 



3 3 



n-n- 1 - 
n— § - 

n-n 



(3.66) 
(3.67) 
(3.68) 
(3.69) 

(3.70) 
(3.71) 

(3.72) 



Note that we can make the weights of the N 5 and Nq seams real by changing the phase of the 
fusion vectors. Consequently, the corresponding partition functions are ambichiral since the 
transfer matrices contain only real entries and thus their eigenvalues must occur in complex 
conjugate pairs. Note also that the weights of the N$ seam can be transformed to be either 
1 or 0. Thus the N 5 seam is seen to be identical to the automorphism seam (|3.35|) which 
implements the Z 2 involution on the Eq graph. 

For D 2 e cases, the ordinary fusion projectors are decomposed according to 



n. 



N x 



1< s < 2£ - 2 



and 



n 2 g-\ = N 2 £-i + N- 



21 



(3.73) 
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and the a- type seams for N 2 e-i and N 2 £ can be obtained by the same process applied to E$. 
In the _D 4 case, the N 3 and N4 a-type seams 

N 3 : TT = 1, = ]Z[ = e2m/3 ^ 3 - 74 ) 



are identical, up to the phase of the fusion vectors, to the Z 3 automorphism seams and yield 
ambichiral partition functions. 

3.6 Construction of a-type seams for Type II theories 

So far we have discussed only the a-type seam of Type I graphs. On a Type II graph G, 
the seam cannot be associated with N a since these graph fusion matrices contain negative 
integers which fail to give a meaningful description in terms of lattice paths. Instead, for a 
Type II graph, an a-type seam must be associated with a node a G H in the parent Type I 
graph rather than to a node a G G. For the A-D-E graphs, only D 2 £ + i and Ej are of Type 
II. For D 2 e + i we do not need a-type seams only the r- and s-type seams. So here we consider 
only the case of G = E 7 with parent graph H = D w . 

The s-type seams of E 7 constructed via ()3.32|) are labelled by s G Ayj. We decompose 
these according to (|3.73|) 



n 



(E 7 D 10 ) 



l<s<8, and n 9 = nf 7 Dw) +n\* 7Ul0 > (3.76) 

where the intertwiners {ni^ 7 Dl0 ^}a€D w of E 7 relative to D w ()3.25|) form a representation of 
the -D10 graph fusion algebra. Explicitly, we have 



(E 7 D10) 



n 9 



\ 



( 



\ 



( 



n, 



(E7D10) 
9 



+ 



+ 



\ 



n 



(E 7 D 10 ) 
10 



The ni E7Dw) 



(3.77) 



fusion vectors for Ej can be obtained by unitary transformation. However, 
as n 9 is complicated, we need to introduce another unitary transformation formula before 
we proceed. Since (13.49(1 is unitary, its inverse is also unitary. By the same argument, 
suppose that U (e, d) is known for some given s, d, e. Then for any b,cEG with N ab c 7^ 
satisfying 

V0 and n sb c ' n 2d d = n sb c 5 c ^ c (3.78) 



n 2( 

is) a 



the N a fusion vector U 



U„ (e, d) 



(6, c) is given by the unitary transformation 



[b, c) 



(b, c) 




(3.79) 
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and can be expressed as 



(b, c) 



b 7 c 



e fi 



(3.80) 



Comparing ()3.80|) with ()3.51|) we note that it does not contain the crossing factor. Depending 
on the condition ()3.48j) or ()3.78|) . the unknown N a fusion vectors can be obtained from known 
ones by either (J3.51|) or (J3.80|) . 

Note, however, that for (6, c) = (3, 3) in n 9 , there does not exist (e, d) which satisfy 
either (j3.48J) or ()3.78|) . Thus, we need to extend the range of the transformation in order to 
make it unitary 



{ *7 (9) 10 (2, 2), 10 (2, 4)} W {U? 10 (3, 3), llf "(3, 3)} 
{ (9)9 (2,2), ? 7 (9)9 (2,4)}^ {? 7 (9)9 (3,3), C /(1,3)} 



(3.81) 



T -v(9)o 



and the transformations to U ' (3, 3) for both a = 9, 10 must be in the nullspace of 17(1, 3). 
Thus 



*7 (9)9 (3,3) 



^3(l/>2+V>4) 




U(9)9( 2 ,2) x 




U(9)9(2,4) X 



(3.82) 



where the scalar factors in ()3.82|) cancel the crossing factors in (|3.51j) when we take the 

* (9) 9 ~ (9) 9 

inverse of 17(1,3) i-> U (2,2) and J7(l,3) i-> U (2,4) respectively. 

There are 58 possible spin configurations for non-zero s-type seam weights 9 ^. We 

constructed all the a-type seams for ng ErDl0 ^ and n^ 7Dl0 ^ symbolically and confirmed that 
(|3.45J) is satisfied for all possible spin configurations at the four corners of the seam weight 
and for all possible bond variables. Note that, by the quantum symmetry, the partition 
functions of x = (a, b) = (9, 1) and (3, 1) are complex conjugates. Thus, the spectra of the 
transfer matrices of the respective ng E7Dl0 ^ and seams must be complex conjugates and 
this is verified numerically. 



4 Transfer Matrices 

Given the fusion hierarchy, we build commuting transfer matrices for different fusion levels 
and boundary conditions: on the torus and on the cylinder, with or without seams. 

4.1 Torus transfer matrices 

The transfer matrix for the (p, g)-fused A-D-E lattice model with an (r, s, r])-seam, on a 
toroidal square lattice is given, in the basis of the cyclic paths in iV steps plus the seam, 
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with bond variables between adjacent spins, by the product of the corresponding Boltzmann 
weights. The entries of the transfer matrix with an (r, s, rf) seam are given by 



bi b2 



u,0 \b,f3) 



-otff- 



r K0 



ajv+i 



Pn+2 



UJV43 






b N /3 N 




( 


In 


Tiv+i 




v a^v oat 








( b]y + 2 


lN+2 


) ^M. 1 


lN+2 


a N+2 , 







«,£ I x 

7i 
a 1 



(4.1) 



where the sum is over all allowed vertical bond variables. The usual periodic boundary 
condition is obtained for (r,s,rj) = (1, 1, 1). The s-type seam can be replaced with a single 
a-type seam or a pair of a and b seams. Recall that the a and 6-type seams derive from 
the two braid limits u — > ±ioo respectively of the s-type seams and are related by complex 
conjugation. Indeed the definition can be generalized to accommodate an arbitrary number 
of seams. Because the seam faces, other than the automorhism seams, are modified bulk 
faces they automatically satisfy the GYBE. They can therefore be moved around freely 
with respect to each other and the bulk faces without effecting the spectrum of the transfer 
matrices. However, in the D 2 e cases when there are several seams, their order in general 
cannot be interchanged because the automorphism seams do not commute with the a-type 
seams. 



4.2 Non-commutativity of seams 

To understand the origin of non-commutativity of seams let us begin by considering D4. The 
graph D4 exhibits an S3 symmetry on the external nodes T = {1,3,4}. This permutation 
group contains non-commuting two- and three-cycles 

S 3 = {();(13),(14),(34);(134),(143)} (4.2) 

In the lattice model seams, these symmetries are realized by the Z2 and Z3 automorphism 
seams (j3.35J) with 77 G S 3 . Notice however that the graph fusion matrix implements the 
permutation (13 4) and N4 implements the permutation (14 3). Accordingly, on the lattice 
these automorphisms are implemented by the a-type seams with a = 3,4. This means that 
the fused seams need only be supplemented by the Z 2 transposition a = (3 4) to generate 
all of S3. Moreover, D4 is the only graph with an automorphism group other than Z2, so we 
can always restrict to automorphisms rj = <j K_1 G Z 2 , K = 1,2. 

The fused seams, together with the bulk face weights, satisfy the generalized Yang- 
Baxter equation and thus they commute with each other and can propagate through the 
bulk face weights along the row 0]. Two transfer matrices differing by a propagation of 
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fused seams differ only by a similarity transformation. Thus the spectrum of the torus 
transfer matrices with a given set of regular fused seams does not depend on the order of 
these seams nor on their positions. On the other hand, the automorphism seams commute 
with the bulk face weights and the r- and s-type fused seams but not, in general, with the 
a-type seams related to N fusions. This is manifest in the -D4 case since the transposition 
(3 4) does not commute with the three-cycles (13 4) and (14 3) implemented by N 3 and N 4 . 

Consider several such seams on the lattice. The non-commutativity of seams shows up 
by the fact that we get different resultant seams when placing given seams in different order. 
For instance, if a = (3 4) 

(143)(34)(134)(34) = (134), (1 4 3)(1 3 4) (3 4) (3 4) = () (4.3) 

and accordingly the product of four (r, a, rj) seams 

W(i,4,i) W {ltl>a) W (h3:1) W {1>lj<7) ~ W(i )3il ) (4.4) 
gives the same spectrum as the single seam Wr\. } z,x) whereas the combination of seams 

W {1A:1) W( lt3>1) W (1A>a) W (1>lj<7) ~ W(i,i,i) (4.5) 

yields the modular invariant partition function. Notice that Wn^Wni^ ~ W^i, 1,0-) ^(1,3,1) 
gives the same spectra even though (13 4) (3 4) 7^ (3 4) (13 4) since the positions of the two 
seams can be interchanged by propagating one of them full cycle around the periodic row. 
Consequently, four or more seams are required to see the effects of non-commutativity in the 
spectra. 

The same phenomenon is exhibited for the whole family D 2 e of Type I models. The 
a-type seams (1,2£ — 1,1), (1,2^,1) related to iV^-i an d N 2 e do not commute with the 
automorphism seam (1,1, a) where a = (2£ — l, 21) is the Z2 transposition, although their 
sum does commute as indicated by the relation N-ii-\o = crN^. 

4.3 Integrable seams on the cylinder 

Although twisted partition functions occur on the torus, it is striking to see that the Ocneanu 
algebra still plays a role on the cylinder. Indeed, Petkova and Zuber |25j give the minimal 
conformal partition functions on the cylinder with a seam (r, x) = (r, a, b, k) and boundary 
conditions (r', a) and (r", b) as 

Z(r',a)\(r,x)\(r",b){q) = ^ ( N r"' N ^~ 2 ^'^ ( n ^x)a Xr"',s (<?) (4.6) 
(r"',s) 

Another remarkable observation 0] is that the twisted partition functions on the torus can 
be written as a bilinear combination of cylinder partition functions, summed over some 
boundary conditions. 

In fact, the Ocneanu graph labelling the seams gives a complete set of boundary con- 
ditions, not only on the torus but also on the cylinder for the continuous conformal field 
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theory as well as for the integrable statistical mechanics model. In the latter context, a 
full set of integrable boundary weights of type (r, a) for Type I theories can be obtained 
by propagating a seam of type (r, a) to the boundary and combining it with the simplest 
boundary condition (r, a) = (1,1) called the vacuum. Similarly, for Type II theories, the 
seams (r, c) with c in the parent graph H can be propagated to the boundary to produce 
an integrable boundary condition by combining it with the vacuum. However, in this case, 
only \G\ of the \H\ values for c produce linearly independent boundary conditions. This is 
in accord with the fact that only \G\ of the \H\ intertwiners n c of G relative to H are 
linearly independent 

~ (4.7) 



71 



(GH) 



E 

aeG 



n 1 iV (G) 



In other words, although G-type boundary conditions are applied on the edge of the cylinder, 
only if -type seams propagate into the bulk. With this caveat, the algebra of defect lines (or 
seams) can be applied on the cylinder as well. 

The vacuum boundary condition corresponds to (r,a) = (1,1). The (l,a) boundary 



weights, for two g-adjacent nodes of G, c and a (i.e. 



71 



7^ 0) are given by 



B 



1 



c 



CM I a 




e<r(l,a)| 



%U^(c,ayUl +1 (c,a) 

Ipa 



ft, 



(4- 



In the Type I case, it is obtained by the action of an a-seam on the vacuum boundary 
condition. The full (r, a) boundary weights are then given by the action of an r-type seam 
on the (1, a)-boundary weight. The double row seams are given by two regular r-seams 
sharing the same extra spectral parameter £, placed one on top of the other, with the same 
spectral parameters as the bulk faces appearing in the double row transfer matrix (see (j4.11}0 

d 



B 



(r,a)\ ^ a 



d S 



b (3 




V-u-(9-i)A,,£) 



-)8 



(l,a) 



(4.9) 



and the left boundary weights are simply equal to the right boundary weights after applying 
crossing symmetry. 

These boundary weights satisfy boundary analogs of the bulk local relations. The Gen- 
eralized Boundary Yang-Baxter Equation or reflection equation is 





(4.10) 



f ' f f 

which is proved using the GYBE (I3.16J) and the abelian property f)3.3j) . We refer to [4|, |2 
for the boundary crossing equation. 
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The double row transfer matrix is given by two rows similar to the one appearing in the 
torus transfer matrix, with spectral parameters u for the bottom one and \i — u — (q — 1)A 
for the top one, where /i is a fixed parameter and q is the vertical fusion level. The boundary 
condition is not cyclic but given by the boundary weights ()4.9|) . 



(a, a | T 



,/;,/ 



(rL,a L )\(r,s,r])\(r R ,a R ) 




-(q-l)\) 



a N 



U 



( »~ U V 



a JV+3 




(4.11] 



The GYBE (|3.16|) and other local relations imply that double row transfer matrices with the 
same boundary conditions and boundary fields commute 



(rL,ai,)\(r,8,Tj)\(rR,a.R)\ > S-k) S; (rL,a.L)\(r,s,r))\(r R ,a R ) 

T Ca^i)|(r,*,»j)|(ra,Oii)( V '^'^^) T lk,a. L )\(r,s, V )\(r R ,a R )( U ' ■ ( L..L2) 



5 Finite-Size Corrections and Numerical Spectra 
5.1 Finite-size corrections and conformal spectra 

A critical A-D-E lattice model with a spectral parameter in the range < u < A gives rise 
to a conformal field theory in the continuum scaling limit, namely, an s£(2) unitary minimal 
model. The properties of the A-D-E lattice model connect to the data of this conformal 
field theory through the finite-size corrections to the eigenvalues of the transfer matrices. 

Consider a periodic row transfer matrix T(u, £) with a seam x of type (r, s, rj) or (r, a, if) 
and N faces excluding the seams. If we write the eigenvalues of this transfer matrix as 

T n (u)=exp(-E n (u)), n = 0,1,2,... (5.1) 

then the finite-size corrections to the energies E n take the form 

E n (u)=Nf(u) + f r (u,0 (5.2) 

27T / C — — — — \ 

+ — (J - — + A n + A n + k n + k n ) sin i? + i(A n - A n + k n - cos^J + o 

where f(u) is the bulk free energy, f r (u,£) is the boundary free energy (independent of s, 
a and 77), c is the central charge, A n and A n are the conformal weights, k n , k n G N label 
descendent levels and the anisotropy angle $ is given by 

d = gu (5.3) 
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where g is the Coxeter number. 

On a finite M x N periodic lattice, the partition function can be written as 



Z^ N = exp(-MNf(u)-Mf r (u,0)Z x (g) 

= TiT{u^) M = Y J T n {u) M = Y J e M-ME n {u)) (5.4) 

n n>0 

where Z x (q) is the conformal partition function and q = exp(27ur) is the modular parameter 
with r = ^exp[i(7r — •&)]. Removing the bulk and boundary contributions to the partition 
function on a torus leads to the twisted partition functions Z x (q) 0] described in Sections 12. 31 
andl2~l 



5.2 Bulk and seam free energies 

In ji^, we showed that the row transfer matrix with an (r, s,^)-seam satisfies the inversion 
identity hierarchy 

TlT\ = sl-M^f-ift I + s^fo Tl (5.5) 
where T q k = T-? (it + kX) are vertically g-fused transfer matrices, 

f"= [<(«)]". ■i^.+q ;>i' (5.6) 

In the thermodynamic limit, the second term in (J5.5[) vanishes and the resulting formula is 
called an Inversion Relation. This equation can be solved, using the structure of zeros and 
poles, first at order N and then at order 1 to find the bulk and seam free energies as we 
explain in this section. 

We calculate the bulk and seam free energies, f(u) and f r (v>,£), or equivalently the 
partition function per face k(u) = exp(— f(u)) and partition function per length K r (u,£) = 
exp(— f r (u, £)). Two A-D-E models sharing the same Coxeter number are related by inter- 
twiners so their bulk and seam free energies are the same. Thus we only need to find the 
free energies for the or ABF models [171 ]. 

The bulk free energy f(u) = — log k,(u) for the ABF models was computed by Baxter 

. . f +oc cosh(7r — 2X)t sinhut sinh(A — u)t , 

k(u) = exp / i '- — — i '- dt. 5.7) 

J_ oa tsinli7rt cosh At 

This integral has a closed form when L is even 

L-2 

sin(w + A) A sin(u + (2fc + l)A) 

K(u) = 



sin A ^ sin(w + 2fcA) 
The partition function per face of the Al model satisfies the crossing symmetry 

k(u) = k(X-u) (5.9) 
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and the inversion relation 



k(u) k(u + A) 



sin(w + A) sin(w — A) 
sin 2 A 



(5.10) 



This solution is the unique solution of the inversion relation, crossing symmetry and height 
reversal symmetry which is analytic and non-zero in the analyticity strip Rew G (0, A). 
Likewise, the seam inversion relation for the order one term gives 



K r (u) K r {u + A) 



sin(u + £ + A) sin(u + £ — (r— 1)A) 



sin A 



The range of validity for the parameter £ is 



-A-|<Re(« + 0<-2 



(5.11) 



(5.12) 



Let q be the RHS of (J5.11j) . It is Analytic and Non-Zero in the strip Re« 6 (0, A). 
Furthermore the derivative q' approaches a Constant when Imw -> ±oo (ANZC). Hence we 
can introduce the Fourier transforms of the logarithmic derivatives 



Hk) := I f'Au)e- ku du 

2m J 



d 

du 



0<Reu<A 
+oo 



fr(u) 



T{k)e ku dk 



(5.13) 



(5.14) 



so that (|5.11J) becomes 

F(k)(l + e hX ) 
The solution by inverse Fourier transforms gives 



2m J q{u) 

0<Re (u)<A 



(5.15) 



w 



(5.16) 



0<Re«j<A 



Integrating with respect to u and taking the w integration along the vertical line Re w = e > 
we obtain in the limit e — > 



fr(u) 



log q{u) 1 



log q(u — iw) 



dw 



2 A J_ 00 sinh (L + l)w 
This integral formula admits a closed form under certain conditions 



(5.17) 



K r [U) 



+ 



7-1 

2 

n 



(u + C - 2k\) 



tari i±i( M + e) ll (M + e - ( 2fc-l)A) 
(w + e + A) 



r odd 



L-r 
2 



tan 



L+l, 



u + ^ + A) 



M + e+(2 fc+i)A) 
II — 7 — 7 , c~\ j \ \ — 5 L + r even. 



k=l 



(u + £ + 2k\) 
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The tangent parts are solutions of the homogeneous functional equation and fix up the zeros 
and poles of the sine parts which are solutions of the functional equation (|5.11|) but have an 
unwanted zero in the analyticity strip. 

The seam free energy for s-type seams are given by the braid limit of the r-type seams 
and are constants. The seam free energies of the a-type seams are the same as the s-type 
seam from which they originate. We remove these seam free energies by the normalization 
of the transfer matrices. Lastly, the seam free energy of an 7] seam is zero. 

5.3 Numerical determination of conformal spectra 

The twisted conformal partition functions are obtained numerically from finite-size spectra. 
Since integrable lattice realizations of the s£(2) A-D-E Wess-Zumino-Witten conformal field 
theories are not known we can only obtain numerically the twisted partition functions of 
the unitary minimal A-D-E models labelled by a pair of graphs (A, G). We use Mathemat- 
ica to construct and diagonalize numerically the finite-size transfer matrices T(w, £) with 
specified seams for different numbers of faces N. For simplicity, we restrict ourselves to the 
isotropic conformal point given by u = | and £ = — 3A/2 . For the first ten or so eigenvalues, 
we extrapolate the conformal corrections to N = oo using a combination of polynomial fits 
in the inverse number of faces and van den Broeck- Schwartz [3| sequence extrapolation. 
The resulting sequences give approximations to the rational exponents that appear in the 
g-expansion of the twisted partition functions in increasing powers of the modular parameter 

9- 

In subsequent subsections, we analyse the numerical data for the D4, D 5 ,Dk cases of 
the A-D-E lattice models. The Al cases for L = 3, 4, 5, 6, 7, 9 were reported in 13]. All the 
numerical results confirm the quantum symmetries and twisted partition functions stated in 
Sections 12 . 31 and 12 . 41 Given the coincidence of the construction labels of our integrable seams 
and the conformal labels of Petkova and Zuber, we expect that our list of integrable seams 
will also exhaust the twisted conformal boundary conditions for the exceptional Eq, E-j and 
E 8 cases as well. These cases, however, are too large to convincingly confirm numerically. 

Consideration of a seam of type (r, a) gives access only to a chiral half of the Ocneanu 
graph. In general, to obtain the complete Ocneanu graph, one needs to consider the com- 
position of two seams (r,a,rj) and (r',b,T]'). Since the A-D-E models are labelled by pairs 
(A,G), with the first member always for type- A, it is sufficient to take r' = r. Also, as we 
have seen, we can restrict the automorphisms to a Z2 subgroup of the full automorphism 
group given by rj = a K . Indeed, the only graph G with an automorphism group larger than 
Z 2 is _D 4 and, in this case, the three-cycles of S 3 are reproduced within the N graph fusions. 
We conclude that it suffices to consider integrable seams of the form x = (r, a, b, cr K_1 ). 

For simplicity, we take the absolute values of the eigenvalues which is equivalent to 
taking the modular parameter q real 



| T flV)| =e -"/-/r 



2tt _£_ 
N 12 




) 



) 



(5.19) 
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We thus numerically estimate the conformal dimensions or exponents 

x n = A + A + k n + k n (5.20) 

and ignore the spins s n = A — A + k n — k n . Since we are at the isotropic point the geometric 
factor sin-i? = 1. To obtain the bulk free energy, we extrapolate the sequence — jt 
The seam free energy is obtained, in a similar way, after removing the bulk contribution in 
(j5.19j) . An s- or a-type seam contributes a constant to the seam free energy which is removed 
by our choice of normalization. The extrapolated numerical values for the bulk and seam 
free energies agree with the analytic results (j5.7|) and (j5.17|) within an accuracy of ±0.3%. 
To estimate the exponents x n , we extrapolate the sequences 

^ = ~ (log \T^\ +Nf + f r - (5.21) 

and compare values and degeneracies with g-series of the twisted partition functions Z x (q) 


The accuracy of our numerical results is restricted by the data for different system 
sizes iV which, in turn, is limited by available computer memory. A typical maximum size 
matrix that we can construct and diagonalize is around 4500 x 4500. The dimension of 
a transfer matrix with N faces and an (r, a) seam is Tr (n^n r Na). This grows rapidly as 
the number of nodes in G increases and grows exponentially with N. In practice, this 
means we are typically restricted to \G\ < 6 and to system sizes N < 12. Furthermore, 
because of parity constraints, we are either restricted to odd or even system sizes N so we 
are typically extrapolating sequences of length six. Nevertheless, because of the quantized 
values of the conformal weights, the integer spacing of conformal towers and recognizable 
degeneracies of the characters, we are able to identify the various twisted partition functions 
with considerable confidence. 

5.4 Numerical spectra of (A4, D4) 

The (A4, D4) lattice model corresponds to the critical 3-state Potts model and is special in 
that the -D4 graph admits an S3 automorphism group. 

However, we only need seams of the form (r, a, k) = (r, a, o^" 1 ) with k — 1,2 and 
a = (34) the Z 2 transposition represented by 

(5.22) 

To illustrate our numerical procedure consider the twisted partition function Z( lj2 ,i){q) 
with seam (r,a,n) = (1, 2, 1). The g-series at the isotropic conformal point with real q is 

Z(i,2,i)(g) = Xi, 2 (Xi,i + 2x li3 ) + X 3 , 2 (X3,i + 2 X3,3) 

= (Xi,2 + Xi,J (Xi,i + Xi,5 + 2 Xi, 3 ) + (X 3 ,2 + Xm) (X 3 ,i + X 3 ,5 + 2 X3,3) ( 5 -23) 

= g-l/15(2g n /120 + g l/8 + g 17/40 + 2 g 71/120 + 2g 19/24 + ? 37/40 + 4g 131/120 + o ^131/120\j 
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Figure 6: Energy levels of the minimal (Ai, D±) model with (1, 2, 1) seam. The numerical 
values of Xn are plotted against 1/N. The polynomial fits intersect the y-axis at the 
extrapolated values x n . Notice that eigenvalues can cross for small values of N. 



so that the exponents, counting degeneracies, are given by the sequence 

Ir \ — f 11 JJ_IiZJLJL191937131131 \ (r oa\ 

X^ni L 120 ' 120 ' 8 ' 40 ' 120 ' 120 ' 24 ' 24 ' 40 ' 120 ' 120 ' ' ' ' J yo.^j 

A plot of our numerical data is shown in Fig. El Notice that some eigenvalues cross so 
that it is not possible to simply order the eigenvalues according to their magnitudes at a 
given value of N and this complicates the extrapolation procedure. Although we have not 
systematically done so, it is possible to remove this ambiguity in the identification of each 
eigenvalue at a given N by examining the patterns of the zeros of the eigenvalue Tn (u) in 
the complex u-plane as explained in [2^| . The theoretical data (|5.24jl is to be compared with 
our numerical values of {x n } for the (1, 2, 1) seam shown in Table [7| The agreement over 
the first 10 levels is certainly good enough to unequivocally identify the (1,2,1) integrable 
seam as giving rise to the (1, 2, 1) twisted partition function. 

Typical numerical data for two other seams is also shown in Table |7| In this way we 
can identify all of the integrable seams with the corresponding twisted boundary conditions 
in Table El given by Petkova and Zuber jjj. In particular, we confirm that our integrable 
seams give a realization of the complete set of conformal twisted boundary conditions for 
the 3-state Potts model. 



5.5 Numerical spectra of (^6, D§) 

The numerical spectra of the minimal (Aq,D§) model can be obtained similarly to the 
(At, Da) model. The graph D 5 , however, is a Type II graph with parent graph A7. The 
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(1,2,1) 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Exact 


11 

120 


11 

120 


l 

8 


17 

40 


71 

120 


71 

120 


19 

24 


19 

24 


37 

40 


131 

120 


Num. 


0.0915 


0.0915 


0.1250 


0.4258 


0.5934 


0.5934 


0.7969 


0.7969 


0.9215 


1.0421 


|diff.| 


io- 4 


IO" 4 


3-10" 5 


8-10~ 4 


0.0018 


0.0018 


0.0053 


0.0053 


0.0035 


0.0496 


(1,3,1) 


Exact 


2 

15 


7 

15 


7 

15 


2 
3 


2 
3 


17 
15 


17 
15 


4 
3 


22 
15 


22 
15 


Num. 


0.1335 


0.4658 


0.4658 


0.6681 


0.6681 


1.1640 


1.1640 


1.3364 


1.4445 


1.4445 


|diff.| 


IO" 4 


8-10" 4 


8-10" 4 


0.0014 


0.0014 


0.0306 


0.0306 


0.0030 


0.0222 


0.0222 


(3,1,2) 


Exact 


i 

20 


3 
20 


3 

20 


11 

20 


n 

20 


13 
20 


13 
20 


21 
20 


21 
20 


21 
20 


Num. 


0.0492 


0.1501 


0.1501 


0.5479 


0.5479 


0.6503 


0.6503 


1.0577 


1.0577 


1.0629 


|diff.| 


8-10" 4 


10~ 4 


IO" 4 


0.0021 


0.0021 


3-10~ 4 


3-10" 4 


0.0077 


0.0077 


0.0129 



Table 7: Numerical exponents x n for the (v4 4 , _D 4 ) minimal model with seams of type 
(r,a, K ) = (1,2,1), (1,3,1) and (3,1,2). 

integrable seams are thus labelled by (r, s) G (A e , A 7 ) where the s-type fusion seams are 
subject to the fused adjacency matrices n s of D 5 . Note that, although n 2 =n 6 , the eigenvalue 
spectra associated with these seams off the isotropic point are complex conjugates of each 
other in agreement with the quantum symmetry of the Oceneau graph. Also, note that 
nj and the corresponding seam (1,7,1) = (1,1,2) implement the Z 2 involution. Thus, the 
labels (r, s) suffice to label all the distinct nodes on the Oceneau graph for D§. In fact, this 
is a general feature for D 2 e+i for which the Z 2 involution is implemented by n^-\. Our data 
for seams of type (r, a, k) = (1, 2, 1), (3, 3, 1) and (1, 1, 2) are shown in Table |H1 

5.6 Numerical spectra of (Ag, Dq) 

For D 6 we have fewer data points than for Z} 4 so the precision of the extrapolated values of 
the exponents is not as good. Nevertheless, the characteristic degeneracies of the exponents 
in the twisted partition functions are faithfully reproduced. Our data for seams of type 
(r, a, k) = (1, 5, 1), (8, 2, 1) and (3, 5, 2) are shown in Table M 

As dictated by symmetry we observe generally, for (A 4 £_ 4 , D21) , that the a-type seams 
(1, 21 — 1,1) and (1, 2£, 1) give the same spectra for the transfer matrices. This is in accord 
with the symmetry of the graph fusion matrices. Indeed, if a implements the Z 2 transposition 
for (A^_ 4 , D 2 i) then 

a N 2e -i = N 2i a (5.25) 
so that N 2 e-i is similar to N 2 i under interchange of the nodes 2£ — 1 and 2£. 
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(1,2,1) 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Exact 


23 
224 


23 

224 


2Y 

224 


5 

32 


83 
224 


1U7 

224 


135 

224 


135 

224 


167 

224 


167 

224 


Num. 


0.102 


0.102 


0.1201 


0.1562 


0.3728 


0.481 


0.6074 


0.6074 


0.7527 


0.7527 


|diff.| 


7-KT 4 


7-HT 4 


5-KT 4 


5-10~ b 


0.0023 


0.0033 


0.0047 


0.0047 


0.0071 


0.00714 


(3,3,1) 


Exact 


i 

28 


i 

28 


i 

14 


y 

112 


9 
112 


15 
112 


i 

7 


i 

7 


A 

14 


25 
112 


Num. 


0.036 


0.036 


0.0722 


0.0813 


0.0813 


0.1354 


0.1459 


0.1459 


0.215 


0.2275 


|diff.| 


2-10" 4 


2-10~ 4 


7-10- 4 


9-10" 4 


9-10~ 4 


0.0015 


0.0031 


0.0031 


0.0007 


0.0043 


(1,1,2) 


Exact 


3 
112 


15 
112 


5 

16 


4 
7 


4 
7 


yy 

112 


115 
112 


115 
112 


127 
112 


127 
112 


Num. 


0.027 


0.1347 


0.3126 


0.569 


0.569 


0.8785 


1.0310 


1.0310 


1.1673 


1.1673 


|diff.| 


2-10" 4 


8-10~ 4 


MO" 4 


0.0026 


0.0026 


0.0054 


0.0042 


0.0042 


0.033 


0.033 



Table 8: Numerical exponents x n of the (A 6 , D 5 ) minimal model with seams of type (r, a, k) = 
(1,2,1), (3,3,1) and (1,1,2). 



(1,5,1) 


n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Exact 


2 

45 


2 

15 


4 
15 


8 
15 


8 

15 


29 

45 


29 

45 


14 

15 


47 
45 


47 

45 


Num. 


0.0447 


0.1345 


0.2688 


0.5267 


0.5267 


0.6438 


0.6438 


0.9625 


1.0256 


1.0256 


|diff.| 


2-10~ 4 


0.0011 


0.0021 


0.0066 


0.0066 


7-10- 4 


7-KT 4 


0.0292 


0.0188 


0.0188 


(8,2,1) 


Exact 


13 
120 


13 
120 


43 
360 


17 
120 


7 

40 


41 
120 


151 
360 


61 
120 


73 
120 


73 
120 


Num. 


0.1081 


0.1081 


0.1193 


0.1408 


0.1749 


0.3461 


0.4214 


0.5061 


0.6252 


0.6252 


|diff.| 


2-10~ 4 


2-10~ 4 


M0~ 4 


8-10~ 4 


8-10" 4 


0.0044 


0.0019 


0.0022 


0.0169 


0.0169 


(3,5,2) 


Exact 


l 

20 


l 

20 


i 

12 


5 

36 


5 

36 


7 

36 


49 
180 


49 
180 


79 
180 


79 
180 


Num. 


0.0503 


0.0503 


0.0854 


0.1433 


0.1433 


0.1941 


0.2741 


0.2741 


0.4276 


0.4276 


|diff.| 


3-10~ 4 


3-10- 4 


0.0020 


0.0044 


0.0044 


0.0004 


0.0018 


0.0018 


-0.0112 


0.0112 



Table 9: Numerical exponents x n of the (A 8 , D 6 ) minimal model with seams of type (r, a, k) = 
(1,5,1), (8,2,1) and (3,5,2). 
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6 Discussion 



In this paper we have constructed integrable seams for each twisted conformal boundary 
condition x = (r, a, b, k) of the A-D-E s£(2) minimal models. Our construction labels are 
identical to the conformal labels of Petkova and Zuber. For Ej, we have not succeeded in 
constructing single integrable seams for the special seams denoted (r,X), even though we 
have constructed integrable seams of type (r, 6, 2, 1) and (r, 4, 2, 1) which suffice to give the 
partition functions, so these seams remain somewhat mysterious. 

The construction of our integrable seams involves a new type of fusion, related to the 
graph fusion matrices N a and labelled by the nodes a of the graph G, rather than the 
usual Young tableaux. We find that the left-and right-chiral halves of the Ocneanu fusion 
algebra are related to the braid limits u — > ±ioo of the seam weights. The quantum sym- 
metries and twisted conformal partition functions have been verified numerically for the 
Al, L = 3,4,5,6,7,8,9 and Dl, L = 4,5,6 cases. In general, the numerics reproduce the 
first 10 exponents to about a 1 or 2% accuracy. Together with the agreement of the exact 
degeneracies, the numerics give strong evidence of the identification of our seams with the 
corresponding twisted conformal boundary conditions. The origin of non-commutativity in 
the Ocneanu graph fusion algebra is traced to the existence of graph automorphisms which 
do not commute with the fusions. 

Lastly, we point out that the parameter £ appearing in the integrable seams is an arbi- 
trary complex parameter. We have fixed its value to ensure conformal boundary conditions 
in the continuum scaling limit. By choosing it to have an imaginary part that scales appropri- 
ately with N (the number of columns in the transfer matrix), it is possible ji^ to perturb the 
boundary away from the conformal fixed point. This will induce renormalization group flows 
between the fixed points representing the various twisted conformal boundary conditions. 
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